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*^ It is shown that perturbative reheating can reach a sufficiently high temperature 

00 with small or negligible inflaton decay rate provided that the inflaton potential 

i^ becomes negative after inflation. In our model, inflaton and dark energy particle 

^ are two independent scalar flelds, and, depending on the mass of the inflaton and 

^ its coupling to matter flelds, there is a possibility that the remaining inflaton after 

^ — ' reheating can become a dark matter candidate. 
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Modern cosmology [1, 2] based on the inflation demands three seemingly different types of 
particles, whose identities are still elusive. Inflaton that drives inflation is one, dark matter 
(DM) is another, and the third one is to explain the nature of dark energy. There have been 
attempts to unify these different fields for the sake of simplicity [3, 4, 5]. Some have tried to 
reduce the number of fields by relating some of them. For example, ref. [6] considers a connection 
between the inflaton and dark energy. In this Letter, we propose a possibility to link the inflaton 
and DM^, while introducing a new mechanism of reheating. 

The model we consider has two scalar fields: one is the infiaton that drives inflation, while 
the other is responsible for later expansion of the universe as dark energy. The total potential 
energy takes the form of 

V{<P,ii) = V^i<P) + V2{v)>0, (1) 

where is the inflaton and rj is the (late time) dark energy particle (DEP). carries no Standard 
Model charges and interacts with fermions only in terms of Yukawa couplings of A/0-?/'-?/', while 
T) does not interact even with fermions. Note that and r] do not interact with each other. 
We assume that Vi(0) becomes negative at some point^, and ¥2(1]) becomes dominant at later 
time to sustain the expansion of the universe. We also assume that 77 remains massless to avoid 
any localization of dark energy, while inflaton could become massive. Since we could actually 
take ¥2(1]) to be constant [3, 5], we choose 

V2{V) = -V^l(0min) (2) 

such that 

I > otherwise. 
This choice can actually avoid flne-tuning because the magnitude of the current cosmological 
constant becomes equivalent to the vacuum fluctuation of the inflaton, which can be small 
enough. 

The ideal potential for the inflaton would be of the form 



-^See [7] for a different proposal. 

^A negative scalar potential is also considered in [8, 9], in which the negative part of the potential is identified 
as big crunch before the big bang. In our case, however, there is no crunch because of ij such that V{(j), ij) > 
always. 
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Figure 1: Inflaton potential for eq.(4). 



Figure 2: Exaggerated inflaton poten- 
tial for eq.(4). 



where a,b,c,d,(f)Q are constants, and shown in Fig.l and Fig. 2. This provides a slow-roll 
potential initially, then changes to negative toward a stable vacuum, and, finally, approaches 
zero asymptotically. The parameter 0o controls the number of e- foldings A'^e. What we further 
need is 0mm should be close to 0o- Note that, in arbitrary units, for 6^1, i.e. assuming b is 
small so that the second term is negligible in computing Ne, then 
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(5) 



To get the needed e-foldings Ng ~ 60, for example, for a = 40, we need 0o — 0.21. To get a 
reasonable 0mm close enough to 0o, we need to adjust c, d properly. The natural unit expression 
can be obtained by assigning a suitable energy scale to 0o, etc., but it is not essential here. 

The universe cools off during inflation, so it needs to reheat to a sufficiently high temperature 
in the radiation-dominated era[l]. (See [10] for a recent review.) It is commonly known that 
the usual perturbative reheating cannot generate a sufficient temperature due to normally 
small coupling constant between the inflaton and matter fields. To overcome this, the idea of 
preheating is introduced in [11, 12], but it requires a coupling between the inflaton and another 
scalar fleld.^ However, if one can achieve enough reheating, the preheating is unnecessary 
and an extra scalar fleld can be avoided. In this Letter, we shall propose a new perturbative 
reheating mechanism that can generate a sufficiently high reheating temperature. 

A sufficiently high reheating temperature requires rapid transfer of the inflaton energy to 
radiation right after inflation. There are two ways of achieving this: one is to take advantage 

■^See [13] for a proposal of preheating with fermionic couphng. However, preheating cannot take place until 
the inflaton becomes massive. 



of the coherent nature of the inflaton (near the new vacuum), another is to take advantage of 
rapid growth of inflaton kinetic energy when the slow- roll potential turns steep [14]. Here we 
shall take the latter path, but it is also related to the former implicitly. 

Consider the equation of state for the inflaton, p^ = w^p^, where 

P0 = |0' + l^i(0), (6a) 

- -102-14(0), (6b) 
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SO that the equation of state parameter w^ can be expressed as 



lj.2 



As the inflaton evolves, w^ is not constant but varies. If Vi(0) > 0, 1 > w,^ > — 1 always. 
However, for p^ > 0, if (and only if) Vi(0) < 0, w^ > 1 is allowed. (If p^ < for Vi(0) < 0, 
w^ < —1.) Note that, if Vi(0) < 0, p^ can be as small as possible. Hence, as the inflaton loses 
kinetic energy, i.e. p<^ — > 0, w<^ ^ 1 can happen. 

The possibility of w^ ^ 1 is the key to successful perturbative reheating in our model. 
Hence, we call it ekpyrotic reheating, borrowing the word from [9, 18]. In general, the inflaton 
can lose its energy by decaying into radiation, scattering with radiation, and/or annihilating 
into radiation. With taking all these processes into account as the source of damping, the 
evolution equation for the inflaton can be generalized for arbitrary w^, also using eq.(7), as 

+ {3H + A^)p^ + J:^pI = 0. (8) 



w^ + 1 

A^ is the prefactor of the linear damping term, which includes contribution from scattering as 
well, and the leading contribution is from usual decay processes of such that 

A,/, = r^ + ■ ■ ■ , (9) 

where the detail of the ellipsis for scattering is not important for our purpose. And 

S* - ^ (10) 

for annihilation processes with some energy scale parameter E^f,. If is massive, E^ = m^. If 
is massless, the average value E^ = p^(t)/n^{t) ^ p<f,{ti)/n^{t,i) oc Tj can be taken. In our 
model, the inflaton decay rate is solely given by 

r.-E^-^(i-44J • (11) 



Note that this decay process is dominated by the largest Yukawa couphng constant and smallest 
fermion mass. Typically, if is massive, decay process can occur because T^ ^ 0. However, if 
is massless, decay cannot happen directly and we have to take into account of other processes. 

We shall first demonstrate the ekpyrotic reheating in the case of A^ ^ 0. Since the 
annihilation term can be ignored for A^ ^ 0, eq.(8) can be solved to yield 

P,{t) = pM {"^f^""'' e-^-"'^''^^''-'^', (12) 

where a^ = a(tj). This shows that p^ decreases more rapidly if w^ > 0, compared to the case 
of w,j) = 0. So, we can already anticipate that this will lead to a higher reheating temperature. 
The energy density of radiation p^ satisfies 

pR + AHpR = -p^ - 3H{1 + w^)p^ - Pn = (w^ + l)A^p^ - Pn, (13a) 

H^ = ^{p^ + P^ + Pr), (13b) 

where we assume p^ ^ p^, p^ so that we can ignore p^ terms in the following. Since a{t) ~ 
^2/(3(w+i))^ where ptot = "U^Ptot for the universe, such that 

a 2 1 
H=- = — --, 14 

a 3(w + l)t' ^ ' 

eq.(13a) with eq.(13b) leads to 

Any consistent solution to this equation requires w = 1/3 as it should be since the universe is 
to be dominated by radiation, and the necessary solution is given by 

PRit) = ^^ - pM I e-(-^+i)^^(-*»). (16) 

The first term leads to the maximum temperature for ti ~ Mp\/ ^ p^{ti)^ which is given by 

T^ax~P0(t.)'/'. (17) 

For the reheating temperature, from eq.(12) we obtain 

1 , . 

t 



{w^ + 1)A<^ 

such that 

Trh ~ Vw0 + VMpiA^, (19) 
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which is higher than the reheating temperature for w^ = 0[1, 10]. Note that t ~ 1/A^ is not 
proper in our model because it makes the reheating period rather too long for small A^. Thus 
eq.(18) should be the correct one. Indeed, this reheating temperature can become sufficiently 
high ii w^ ^ 1 even for small A^, hence it justifies the word "ekpyrotic." 

If A^ = 0, the annihilation term can no longer be ignored, then the solution to eq.(8) is 
given by 



P4>{t) P4>{U)\U) ^ 2a -I \u) ' ^^°^ 

where a = {w^ + l)/(u' + 1) and T,^ is taken as the average value. As w^ ^ 1, i.e. a ^ 1, p^(t) 
decreases because t/tj > 1. Using eqs.(13a)(13b) for a :^ 1, in this case, we obtain 

1 1 ^t.V" PSi) (21b) 



QnGiw + lff^ \tj l + |(w + l)S^ti 

For ti ~ M^\j ^ p^{ti)^ the maximum temperature is eq.(17) as before. With 

1 

- ~ S^p^(ti) ~ n4,{ti){av), (22) 

the reheating temperature is 

Trh ~ (^Pi^j ~ sjM^,n^{ti){av), (23) 

which can be sufficiently high since n^{ti){(7v) is large, particularly, for massless inflaton. Note 
that this is true only if w^ :^ 1, hence ekpyrotic. 

The reheating temperature can be constrained by CMB[15, 16, 17]. Ref.[17] claims Trh ~ 
10^ GeV. In r^ 7^ case, this can be achieved with Yukawa coupling constant A/ ~ 10"'^ for 
massive inflaton with m<^ ~ 1 TeV (see eq.(ll)) even for w^ = 0. For smaller Yukawa couplings, 
sufficiently large w^ can achieve the desired reheating temperature even for massless inflaton. 
This indicates that the reheating temperature given by eq.(19) (or eq.(23)) is actually sufficient 
and preheating is unnecessary in our model. If T^ > Hq ^ 10~^® sec^^ such that the inflaton 
decays entirely during reheating, then cannot be DM. However, for w,^ ^ 1, with sufficiently 
small Yukawa coupling constants, we can achieve F^ < Hq so that there is a possibility that 
the remaining inflaton can become DM after ekpyrotic reheating. 

Let us now consider the case in which the remaining inflaton can become DM after successful 
reheating. Once becomes massive, no longer behaves like a (uniform) fluid such that w^ = 

5 



is possible. In our model, the inflaton DM is basically Yukawa interacting scalar DM. (See [19] 
for another example.) We assume that Yukawa couplings are flavor conserving so that we do 
not have to worry about / — )■ /0 process for massive 0. Then, the relic density can be computed 
in terms of the usual Boltzmann eq. for annihilation/creation process [20] 

h^ + 3Hn^ = -{av){nl-nlQ) . (24) 

Let Y = n^/sE and x = m^/T, where se is the total entropy density of the universe and T is 
the photon temperature, then, generalizing ref. [21] for arbitrary w, the Boltzmann eq. can be 
expressed as 



where 



and 



1/2 hs f, , ldlnhcs\ ,^„ . 
y/g^ V 3 dlnT y 

P = ^cfr(T)^T^ (26b) 

SE = hsiT){l+w)'^T^, (26c) 



, 0, for radiation-dominated era, , , 

- 3 = <^ ' . ' (27) 



w + 1 11, for matter-dominated era. 

Note that the prefactor (1 + w) in s^; is due to se oc {p + p)- 

A relativistic derivation of (av) is given in [21], which reads in the lab frame 

2x f^ 
{^v) = -r72r^ de(TVi,^V~<l + 2e)K^{2xVTT~e), (28) 

-"■2l^J Jo 

where e = ^^ — 1 is the average kinetic energy per mass for the Mandelstam variable s, and 
Kn{x) is the modifled Bessel function of the second kind of order n. Note that this formula 
works as long as two incident particles are collinear. If is sufficiently massive as in our case, 
we can use non-relativistic expansion of af lab in powers of e as 

^^lab = E -/' (29) 

n=0 



such that 

(mj) = an 4- 

2 x 



-> = ««+5T + 0(x-)-E^- (30) 



n=0 
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For Yukawa coupling Xf(j)ff, the cross section for 00 — > // in the CE frame is given by 
_ \j Pf 1 



H^ff 16^ p^ Ef 



^ ^ {ml - Am)Y _ 2EJ - (m^ - Am)) ^^ / 2^^m|+2MV 

ml — Am'l'm'j + Am'jEf '^PiPf \ 2-Ef — ml — 2pjj9j 

(31) 
where pf = (Ef — m))^!'^ . Expanding this according to eq.(29), we can obtain 

Q-o — acE,o — -^ 2 ' y^'^i 

Stt m^ 

where «/ = mjlm^. As «/ — )■ 0, ao — )■ and the next nonvanishing expansion coefficient is 

A) 1 

Oi = acE,i — ocE,o = ;; 5"- (33) 

37rm^ 

If is massless, which is relevant for the reheating process, this cross section develops an 
infrared singularity, which can be regularized by introducing a minimal cutoff initial energy. 

With eq.(30), eq.(25) now reads, in the leading order of x, 

<iy ,j .i /""^^'-^ . (34) 

where Feq is to be suppressed by e~^, n accounts the leading nonvanishing term, and 

1 rv 



A. 



n 



:J ^9T^m<p{(^v)n. (35) 
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Taking average value of g^! for an approximation and integrating from xp to Xq, i.e. from the 
freeze-out to the present, we can obtain a solution to this equation as 

™n+l /^"+1 A \ /xn\^+"+''" 

^ K"^ ^ rr^) (I) • <^«' 

where we assume Xq/xf = Tp/T^ ^ 1. 

If the density evolution overlaps from radiation-dominated to matter-dominated era, we 
should solve eq.(34) consecutively from xp -^ Xc ^ xq, where Xc is for the cross-over from 
radiation-dominated to matter-dominated era. Then, the current relic density can be better 
approximated as 



r-^ 



Yq \ Yf n + lj \XpJ \x, 

where Xq/xc = Tc/Tq ~ 10"^. 



If {av) is large enough, the freeze-out temperature can be determined by 



H{tF) = nEQ{(rv), 



where, for large x, 



«EQ 



9fT' 



x^Koix) 



9fT' 



x^l\-\ 



(3J 



(39) 



27r2 ^^ ' 27r2 V 2^ 

Then, assuming freeze-out takes place during the radiation-dominated era, the freeze-out value 
Yp is given by 



^F — ^eqI^^f) 



n^qixp) 15 x/g^^^^iTp) 



X 



l+n 



SEixp) \/27r KsiTp) Mpim^{av)r 



(40) 



From eq.(36) and eq.(40), assuming the density has evolved as if the entire period were 
radiation-dominated era, we can obtain 



Yn = B- 



X 



l+n 



where 



B 



Mpim^{av)n' 



15 



QcsiTp 



(41) 



(42) 



If we use eq.(37) to compensate the density evolution during the matter-dominated era. 



Yn = B 



Xr 



X 



l+n 



(43) 



xoj Mpim^{av)n' 

Note that Yq in this approximation is 10"*^ times Yq of approximation assuming only radiation- 
domination. Then, with Ks{To) = 3.909, Tq = 2.725 K, and Ho = 2.133 x 10~^^h GeV, (^dm is 
given by 

f^DM/i' = 2.051 X 10^ -^Fo- (44) 



Comparing this to the current measured value I^dm^^ = 0.1126 ± 0.0036[22], we can obtain 
conditions on m^ and A/. We will come back to this for a specific case later. 

The above works only if the annihilation cross section is large enough to satisfy eq.(38). If it 
is too small so that eq.(38) cannot be satisfied at all, n-EQixp) should be used to determine the 
relic density. The freeze-out parameter xp for (av) ^ can be determined from the following 
equation [21]: 



Kox) 



dlnT 



For large x, up to x , 



^i(a;) A 3 \ / 15 105 \ ^ ^ 3 345 



7^2 (x) V SxJ \ 8x 128xV 2x 128x 

then eq.(45) leads to an equation for xp 

9 3 345 ^ 

2^ + 128 ~ ^"""^ " ' ^ ^ 

where 

_ dlnhesjT) 

^ = dT ^^^^ 

can be approximated by taking a suitable numerical value. With eqs.(36)(39) and Cy, = 0, 

«0M.^ - ^™.«.(« ^ 7,814 X 10-^^4'V■^ (49) 

where Tq = 2.35 x lO^^^ GeV. Using hsiTo) = 3.909 and /icfflTp) ~ 10, ^ = 0.83 GeV~^ and 
xp — 30 for m^ ^ 1 TeV. 

Note that — > // requires m^ > 2mf, i.e. 1/2 > af, while (fxp — )> // requires m<^ > rrif i.e. 
1 > af. So, there are three possibilities: 

af>l: a = 0, r^ = 0, (50a) 

1 > a/ > I : CT ^ 0, r^ = 0, (50b) 



l>af: a ^ 0, r^ ^ 0. (50c) 

Since we expect m^ > 1 TeV, «/ > 1 for all flavors cannot happen. So, the fate of the inflaton 
depends on these mass ranges as well as Yukawa coupling constants and there are three different 
cases: 

• Case I: The inflaton is not a DM candidate and decays entirely before present day. 

• Case II: The inflaton is a DM candidate and the relic density is given by eq.(49). 

• Case III: The inflaton is a DM candidate and the relic density is given by eq.(44). 

If there is at least one light fermion flavor with Aj > 10^^^ (i.e. F^ > 10^^^ sec^^), it is 
Case I. If \f < 10^^® for all flavors, it is Case II. Case III is the case with some flavors with 
1 > a/ > ^ and other flavors with | > a/. A good example of Case III is the technicolor (TC) 
model with Yukawa interacting sterile scalars[23], in which technifermions have Atc > 10~^^ 
while quarks and leptons have Xq^l ^ 10~^^. Another possibility of Case III without TC is to 
allow the fourth generation of suitable mass with A4 > 10~^^. 
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Let us look into the TC case more in detail. For m^ ~ 1 TeV, at least the lightest 
technifermion mass should be between 1 TeV and 500 GeV. Then, the density evolution is 
dominated by the annihilation process so that we can use eq.(43) to determine the relic density. 
For example, if we assume there is only one technifermion within the mass range with mass 
about 600 GeV, then A/ ^ 3.4 x 10"^ leads to m,^ ^ 1 TeV with xp — 18. If we use eq.(41) just 
for radiation-dominated era, we get A/ ~ 0.38 with xp — 28. Since the difference is significant 
enough, we believe that the relic density is better approximated by eq.(43). This is because a 
significant part of relic density evolution takes place during the matter-dominated era so that it 
should not be ignored. If we have more technifermions within the mass range, smaller Yukawa 
coupling constants can lead to m^ ~ 1 TeV. 

In this Letter, we have demonstrated that a successful cosmological model can be achieved 
with just two sterile scalar fields. In cases II and III, the remainder of the infiaton after reheating 
can become a DM candidate. In this context, it may be possible that massive infiaton can 
trigger local clustering of matter even during the radiation-dominated era, hence it can become 
an early seed of structure formation. The other sterile scalar field as DEP behaves more like a 
free particle and is uniformly distributed with its fiat potential, whose magnitude is the same 
as the absolute value of the infiaton potential energy at a new vacuum after infiation, so that 
the infiaton's vacuum fiuctuation after reheating can account for the magnitude of dark energy. 
We leave the computation of the amount of this vacuum fiuctuation as a future task. 

The infiaton potential we have used has an analogous form to the Fermi surface so that it 
could indicate a more fundamental origin of the potential. Also, even though we have not used 
the oscillatory nature of the infiaton field exphcitly (eq.(8) implicitly uses it[l]), we suspect 
there might be a connection between our ekpyrotic reheating and the fermionic preheating 
considered in ref.[13] to a certain level. It will be interesting to see if this is indeed the case. 

I appreciate Bob Scherrer for many informative conversations, reading the manuscript and 
helpful comments. 
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